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CURVE INTERPOLATION METHOD 
BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to a curve interpolation method for 
obtaining a smooth curve based on data of a sequence of command points and 
performing interpolation on the smooth curve in machining a workpiece along 
a curved surface using a numerical controller. 

2. Description of Related Art 

In machining a workpiece along a curved surface of an object such as a 
mold by a numerical controlled machine tool, etc. according to data of a 
sequence of command points which are created by a CAD/CAM device or a 
profiler device, a curve interpolation is performed based on the data of the 
sequence of command points. For example, there is disclosed a method of 
creating an approximate spline curve based on the data of the sequence of 
command points and perfonning the curve interpolation on the created 
approximate spline curve in JP 2-1 13 305 A. 

In JP 10-240328 A, there is described an interpolation method in which 
vectors of line segments between adjacent twos of command points are 
obtained and a modification amount for each command point is obtained such 
that the sum of squares of differential vectors between the obtained vectors of 
the line segments is the least so as to modify each command point. In this 
document, it is described to set an evaluation range in a sequence of the 
command points and create an approximate curve for successive points in the 
evaluation range, so that the command points are modified based on 
differences between the approximate curve and the respective command points 
in the evaluation range. 

* 

In the CAD/CAM device, a tolerance is set for a target curve created by 
the CAD device and a plurality of line segments are created by the CAM 
device within the tolerance to obtain data of end points of the respective line 
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f(t) x = A x t 3 + B x t 2 + Cxt + D x 

f(t)y = Ayt 3 + Byt 2 + Cyt + Dy 

fCOz^Az^ + B^ + Czt + D; 

Thus obtained curve Ce is shown in FIG. 20. 

Then, interpolation is performed on the defined curve Ce with a unit 
not greater than a set unit in preparing the sequence of command points (Step 
S16). 

In the foregoing embodiment, at the start of the procedure, all of the 
command points P0, PI, P2, ... , Pn-1, Pn are read at Step SI. Alternatively, 
only the necessary command points may be read and the procedure may be 
carried out on the read points, so that the approximate curve is successively 
created while reading the data of the command points to expedite the 
procedure. 

In obtaining interpolation points, i.e. shape-defining points, respective 
two points are interpolated between adjacent twos of the command points in 
the foregoing embodiment, respective points more than two may be 
interpolated between adjacent twos of the command points. Further, in 
creating the approximate curve Cm, the shape-defining points not greater than 
two are selected before and after the shape-defining point Qi. The 
shape-defining points greater than two may be selected. Furthermore, one or 
more of the command points P0, PI, P2, ... , Pn-1, Pn may be used as the 
shape-defining points with the interpolation points Ql, Q2, ... , Q2n. 

If a line segment connecting adjacent two of the command points P0, 
PI, P2, ... , Pn-1, Pn is shorter than a reference value, an interpolated point Pj' 
such as a middle point between such adjacent command points Pj, Pj+1 may be 
regarded as a substitute command point for the adjacent command points Pj, 
Pj+1 which are to be deleted. In this case, it is determined whether or not a 
distance between the adjacent command points Pi and Pi+1 is not greater than 
the predetermined at Step S2, and if the distance is not greater than the 
predetermined value, the above procedure is performed to define a substitute 
command point for obtaining the interpolation points. 
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segments as a sequence of command points to be outputted to the numerical 
controller. 

As shown in FIG. 1, line segments LO, LI, L2, ... are created by the 
CAM device within a tolerance width "2w" set by a tolerance amount "w" on 
both sides of an original target curve Cs which are created by the CAD device, 
and data of points PO, PI, P2, ... at both ends of respective line segment LO, LI, 
L2, ... are outputted to a numerical controller as data of a sequence of 
command points. Since a curve Ce for the curve interpolation is defined 
based on the position data of the sequence of command points, the curve Ce 
may exceed the tolerance width 2w set to the original curve Cs. 

According to the method disclosed in JP 10-240328A, the positions of 
the command points are modified but there is a possibility of creating the curve 
Ce same as that in FIG. 1 based on the modified command points, and it is not 
assured that the curve Ce is created within the tolerance width 2w set for the 
original curve Cs. 

It is probable that the sequence of command points are positioned near 
ends of a band of the tolerance width 2w set to the original curve Cs. 
Therefore, if the curve Ce is defined based solely on the sequence of command 
points, the curve Ce may be positioned considerably remote from the original 
curve Cs. For instance, in the case where the original curve Cs is a circular 
arc, as shown in FIG. 2, the command points PO, PI, P2, ... defined by the line 
segments LO, LI, L2, ... within the tolerance width 2w are positioned remote 
from the original curve Cs by an mount approximately equal to the tolerance 
amount w. In this case, two points Ql and Q2 interpolated on each of the line 
segments LO, LI, L2, ... at a ratio of approximately 0.15: 0.7 : 0.15 are 
positioned on the original curve Cs, as shown in FIG. 3. 

The above circumstance is the same in the case where the positions of 
the command points are modified. For example, in the case where the 
original curve is a circular arc, since a sequence of command points are aligned 
along a circular arc, if a curve approximating the command points is created, 
the sequence of command points are not substantially modified. Thus, the 
modified sequence of command points are positioned remote form the original 
curve Cs by the tolerance amount w. 
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Since the curve Ce is defined to pass the sequence of command points 
PO, PI, P2 if there is an error in calculation for obtaining the sequence of 
command points PO, PI, P2, ... or in approximation by a set unit in the 
CAD/CAM device, the error influences the definition of the curve Ce to lower 
machining precision and cause a vibration of the machine. Thus, the 
definition of the curve Ce to pass the sequence of command points may cause 
deterioration of the precision of the machined surface. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide a curve interpolation 
method capable of obtaining a curve approximating an original target curve 
within a tolerance set for the original target curve based on a sequence of 
command points, and performing interpolation on the obtained curve. 

As describe, the command points PO, PI, P2, ... are positioned within 
the tolerance width 2w and the line segments LO, LI, L2, ; v connecting the 
adjacent twos of the command points (PO, PI), (PI, P2), ... are positioned 
within the tolerance width 2w, as shown in FIG. 5. Therefore, points Ql, Q2, 
Q3, ... interpolated on the respective line segments LO, LI, L2, ... are 
positioned within the tolerance width 2w. According to the present invention, 
a smooth curve Ce approximating an original curve Cs is defined using the 
interpolated points Ql, Q2, Q3, and a curve interpolation for machining a 
workpiece is performed on the smooth curve Ce by a procedure comprising the 
following steps (l)-(7). 

(1) Interpolation points Ql, Q2, Q3, ... are defined between adjacent 
twos of the command points (PO, PI), (PI, P2), ... , as shape-defining points. 
The shape-defining points Ql, Q2, Q3, ... are positioned inner the band of the 
tolerance width 2w, i.e. closer to the original curve Cs than the command 
points PO, PI, P2, ... . 

(2) One shape-defining point Qi and adjacent shape-defining points in 
front of and in the rear of the one shape-defining point Qi to be surrounded by 
these adjacent points are selected. 
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(3) A curve Cm approximating the selected shape-defining points is 
created, as shown in FIG. 6. 

(4) The one shape-defining point Qi is moved towards the 
approximate curve Cm to define the modified command point Qi\ 

(5) Repeatedly executing the above steps (2)-(4) successively on the 
other shape-defining points to obtain the modified shape-defining points Ql f , 
Q2',Q3', ... . 

(6) A curve Ce passing the modified command points Q0 ! , Ql f , Q2 f , ... 
is defined. 

(7) The interpolation is performed on the obtained curve Ce. 

In the step (1), one or more of the command points PO, PI, P2, ... may 
be used as the shape-defining points in addition to the interpolated points QI, 
Q2, Q3, .... In this step (1), the interpolation points Ql, Q2, Q3, ... are 
created with a unit not greater than a set unit in preparing the sequence of 
command points PO, PI, P2, ... 

In the step (1), if a distance between adjacent two of the command 
points Pj, Pj+1 is shorter than a first reference value, an interpolation point Pj 1 
between the adjacent two command points Pj, Pj+1 is used as a substitute 
therefor, as shown in FIG. 7. The interpolated point may be a middle point of 
the adjacent two command points Pj, Pj+1 . 

In the step (3), the approximate curve Cm may be created to have the 
least sum of squares of distances from the selected shape-defining points. 

In the step (4), an amount of moving the one shape-defining point Qi 
for modification may be restricted to a set value. 

In the step (4), if a distance between adjacent two of the modified 
shape-defining points Of, Qj+1' is shorter than a second reference value, an 
interpolation point Qj" between the adjacent two modified shape-defining 
points Qj 1 , Qj+l f is used as a substitute for the two modified shape-defining 
points, as shown in FIG. 8. The interpolated point Qj" may be a middle point 
of the adjacent two modified shape-defining points Qj f , Qj+r. 

In the step (4), the one shape defining point Qi may be moved gradually 
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with a unit not greater than a set unit in preparing the sequence of command 
points PO, P1,P2, .... 

In the step (6), first-order differentiate values Q0 (1 \ Ql (1 \ Q2 (1)l , ... of 
the approximate curve Cm at points corresponding to the modified 
shape-defining points QO', Ql f , Q2 ? , ... are used in creating the curve Ce 
passing a sequence of modified shape-defining points Q0 f , Qf , Q2 f ? . . as 
shown in FIG. 9. 

In the step (6), the curve Ce passing the modified command points Q0' 5 
Ql', Q2\ ... may be a NURBS curve or a spline curve. 

In the step (7), the interpolation on the obtained curve Ce is performed 
with a unit not greater than a set unit in preparing the sequence of command 
points P0,P1,P2, .... 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram showing a conventional interpolation 
method in which a curve is created based on a sequence of command points; 

FIG. 2 is a schematic diagram showing an original target curve and a 
created curve in a case where the original target curve is a circular arc; 

FIG. 3 is a schematic diagram showing intersection between the 
original target curve and line segments connecting adjacent twos of the 
command points in the conventional interpolation method; 

FIG. 4 is a schematic diagram showing an example of an inappropriate 
curve created from the command points according to the conventional method; 

FIG. 5 is a schematic diagram showing interpolation points used as 
shape-defining points according to the present invention; 

FIG. 6 is a schematic diagram showing movement of the 
shape-defining point toward the approximate curve to obtain a modified 
shape-defining point; 

FIG. 7 is a schematic diagram showing substitution of command points; 
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FIG. 8 is a schematic diagram showing substitution of modified 
shape-defining points; 

FIG. 9 is a schematic diagram showing a first-order differentiated value 
for a modified shape-defining point; 

FIG. 10 is a block diagram of a numerical controller for carrying out 
the method of the present invention; 

FIG. 1 1 is a schematic diagram showing a part of a machining program 
for executing the curve interpolation according to the present invention; 

FIG. 12 is a schematic diagram of a sequence of command points for 
which automatic determination of carrying out the curve interpolation of the 
present invention is perfomied; 

FIG. 13 is a schematic diagram of a sequence of command points for 
which another automatic determination of carrying out the curve interpolation 
of the present invention is performed; 

FIG. 14 is a part of flowchart of processing for the curve interpolation 
according to the present invention; 

FIG. 15 is a continuation of the flowchart of FIG. 14; 

FIG. 16 is a schematic diagram of a sequence of command points; 

FIG. 17 is a schematic diagram showing interpolation points between 
adjacent twos of the command points for obtaining shape-defining points; 

FIG. 1 8 is a schematic diagram showing an approximate curve for the 
shape-defining points for obtaining a modified shape-defining point and a 
first-order differentiated value therefor; 

FIG. 1 9 is a schematic diagram showing a substitute point to be 
substituted for adjacent two modified shape-defining points and a substitute 
first-order differentiated value therefor; and 

FIG. 20 is a schematic diagram showing an example of a smooth curve 
created according to the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
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FIG. 10 shows a numerical controller 100 for carrying out the curve 
interpolating method of the present invention. 

A processor (CPU) 1 1 for generally controlling the numerical controller 
100 is connected with a memory 12 including a ROM; a RAM and a CMOS 
memory backed up by a battery; an interface 13 to which a data input/output 
device 34 is connected; an interface 16 to which a display/MDI unit 30 is 
connected; an interface 17 to which an operation panel 31 is connected; a PC 
(programmable controller) 14; a display/MDI unit 30; axis control circuits 18; 
and a spindle control circuit 20, through a bus 22. 

The axis control circuits 18 output commands to the servo amplifiers 19 
for respective axes in accordance with motion commands from the processor 
1 1 . The servo amplifiers 19 drive servomotor 32 for the respective axes in 
accordance with the commands from the axis control circuits 18. Each of the 
servomotors 32 has a position/velocity detector (not shown) for feeding back 
position/velocity feedback signals to the associated axis controller 1 8 for 
feedback control of the position/velocity of the axis. 

The spindle control circuit 20 outputs a spindle speed signal to the 
spindle amplifier 21 based on a spindle speed command from the processor 1 1 
and a speed feedback signal from a position coder provided at the spindle 
motor 33. The spindle amplifier 21 drives the spindle at the commanded 
speed in accordance with the spindle speed command. 

A machining program including data of the command points created by 
the CAD/CAM device or a profile device is inputted into the numerical 
controller 100 from the data input/output device 34 through the interface 13 
and stored in the nonvolatile section of the memory 12. 

G-codes commanding a start of defining a smooth curve and 
interpolation thereon, and commanding discontinuation of the command are 
prepared. A G-code "G06.3" for commanding a start of defining a smooth 
curve and interpolation thereon is programmed at the beginning of data of a 
sequence of command points, and a G-code "G06.9" is programmed at a 
position where this command is to be discontinued in the machining program. 
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The procedure of defining a smooth curve and interpolation thereon is carried 
out in a section between the codes "G06.3" and "G06.9". 

The section in which definition of a smooth curve and interpolation 
thereon are performed may be automatically determined. In this case, the 
processor 1 1 reads the sequence of command points from the machining 
program and determines whether or not the creation and interpolation of the 
smooth curve is to be performed based on angles between adjacent ones of the 
line segments or lengths of the line segments interconnecting the adjacent ones 
of the command points. As shown in FIG. 12, angles between adjacent line 
segments interconnecting adjacent ones of command points PI, P2, P3 and 
command points P6, P7, P8 are smaller than a reference value. Thus, it is 
determined that a smooth curve is to be created and interpolated in these 
sections. Contrary, an angle a between a line segment interconnecting 
command points P3 and P4 and a line segment interconnecting command 
points P4 and P5 are relatively large. Similarly, an angle p between the line 
segment interconnecting command points P4 and P5 and a line segment 
interconnecting command points P5 and P6 is larger than the reference value. 
Thus, it is determined that a smooth curve is not to be created and interpolated 
in the section between the command points P4 and P5. Namely, if an angle 
between the line segments interconnecting the adjacent command points is not 
greater than the reference value, it is automatically determined that a smooth 
curve is to created and interpolated in the section of the line segment, and if the 
angle between the line segments interconnecting the adjacent command points 
is less than the reference value, it is automatically determined that a smooth 
curve is not to created and interpolated. 

Further, as shown in FIG. 13, distances between the adjacent command 
points from the command point PI to the command point P4, and from the 
command point P5 to the command point P8 are relatively short, and a distance 
between the command point P4 and the command point P5 is relatively long. 
Thus, if the distance between the adjacent command points is not greater than a 
reference value, it is determined that a smooth curve is to be created and 
interpolated in the section defined by such command points, and if the distance 
between the adjacent command points is greater than a reference value, it is 
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determined that a smooth curve is not to be created and interpolated in the 
sections between such command points. 

The processing for defining a smooth curve and interpolation thereon 
will be described referring to flowcharts of FIGS. 14 and 15. 

The processor 1 1 reads a sequence of command points PO, PI, Pn-1, 
Pn from a machining program (Step SI). An example of the read command 
points PO, PI, Pn-1, Pn are shown in FIG. 16. 

Respective two interpolation points (Ql, Q2), (Q3, Q4), ... , (Q2n-1, 
Q2n) are defined between the adjacent two command points (PO, PI), (PI, 
P2), ... , (Pn-1, Pn) (Step S2). In this example, an interpolation ratio is set to 
1: 0.7 : 0.15, as shown in FIG. 17. In particular, the interpolation ratio of 
interpolation points (Q2i+1, Q2i+2) on a line segment (Pi, Pi+1) is set such that 
a length of a line segment (Pi, Pi+1): a length of a line segment (Q2i+1, 
Q2i+2): a length of a line segment (Pi, Q2i+1) = 1: 0.7 : 0.15. 

A first shape-defining point Q0 is set to PO and a last shape-defining 
point Q2n+1 is set to Pn (Step S3). Thus, a sequence of shape-defining points 
comprising a sequence of interpolation points Q0, Ql, Q2, ... , Q2n, Q2n+1 are 
created. 

Then, the index i is initially set to "1" in Step S4. A value of the index 
i is determined in Steps S5 and S6, and if it is determined that the index i is f T 
the shape-defining points Qi-1, Qi, Qi+1, Qi+2 (=Q0, Ql, Q2, Q3) are selected 
(Step S17), if it is determined that the index i is "2, 3, or 2n-l M , the 
shape-defining points Qi-2, Qi-1, Qi, Qi+1, Qi+2 are selected (Step S7), and if 
the index i is determined as ,f 2n M , the shape-defining points Qi-2, Qi-1, Qi, 
Qi+1 (=Q2n-2, Q2n-l,Q2n, Q2n+1) are selected (Step SI 8). If the index i is 
determined as "2n+l M , the procedure as described later is performed. 

An approximate curve Cm is created based on thus selected 
shape-defining points by the least-squares method. Namely, the curve Cm 
having the least sum of squares of distances from the selected shape-defining 
points is defined (Step S8). An example of the approximate curve Cm is 
shown as the one-dotted chain line in FIG. 18. The shape-defining point Qi 
designated by the present value of the index i is moved towards the 



approximate curve Cm gradually with a unit not greater than a set unit in 
preparing the sequence of command points and within the tolerance width 2w 
to define a modified shape-defining point Qi 1 with respect to the shape-defining 
point Qi (Step S9). A first-order differentiated value Qi (1)? of the modified 
shape-defining point Qi ! on the approximate curve Cm is obtained and stored 
(StepS 10). 

The index i is incrementally increased by "1" (Step SI 1), and it is 
determined whether or not the value of the index i exceeds lf 2n+l " which is the 
index of the last shape-defining point Q2n+1 (Step SI 2). If it is determined 
that the index i does not exceeds "2n+l M , the procedure returns to Step S5 to 
repeat the processing of Step S5 and subsequent Steps. 

When it is determined that the index i reaches "2n+l M in Step S5, a 
modified shape-defining point Q2n+l f is set as the command point Pn, and an 
orientation vector from the modified shape-defining point Q2n f to the 
command point Pn is stored as a first-order differentiated value Q2n+l (1)f with 
respect to the modified shape-defining point Q2n+r (Step SI 9). Then, the 
index i is incrementally increased by "1" (Step SI 1). As a result, it is 
determined that the index i exceeds "2n+l" which is the index of the last 
shape-defining points Q 2n+l (Step SI 2), and the procedure proceeds to Step 
S13 where the modified shape-defining point Q0' is set as the command point 
P0, and an orientation vector from the command point P0 to the modified 
shape-defining point Ql f is set to a first-order differentiated value Q0 (1)l with 
respect to the modified shape-defining point Q0 f (=P0), With the above 
processing, the modified shape-defining points Q0 1 , Ql f , Q2\ ... , Q2n f , Q2n+l f 
and the first-order differentiated values Q0 (1> , Ql (1 \ Q2 (1 \ Q2n (1 \ 
Q2n+l (1) ' with respect to the shape-defining points Q0, QI, Q2, ... , Q2n, 
Q2n+1 are calculated. 

Further, in this embodiment, as shown in FIG. 19, if there are two 
adjacent modified shape-defining points Qj f and Qj+T having a distance 
therebetween shorter than a reference value in the modified shape-defining 
points Q0 f , QI 1 , Q2n+l f , an interpolation point is created between the two 
points Qj' and Qj+l f . In this embodiment, a middle point of the two points Qj ? 
and Qj+r is used as a substitute for the modified shape-defining points Qj f and 
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Qj+1' which are to be deleted. Also, an average value of the first-order 
differentiated values Qj (1) ' and Qj+l (1)f for the modified shape-defining points 
Qj f and Qj+1' is set as a substituted first-order differentiated value Qj (1) ' for the 
substituted modified shape-defining point Qj', and the original first-order 
differentiated values Qj )f and Qj+l (1) ' are deleted. The first modified 
shape-defining point QO' and the last modified shape-defining point Q2n+1' of 
the sequence of modified shape-defining points and the first-order 
differentiated values Q0 (1)f and Q2n+l (1) ' for the first and last modified 
shape-defining points Q0 ! and Q2n+l f are not deleted. 

With the above procedure, the modified shape-defining points QO 1 , 
Ql 1 , ... , Q2n+1 ? and the first-order differentiated values Q0 (1 \ Ql (1 \ ... , 
Q2n+l (1)f for the modified shape-defining points QO', Ql 1 , ... , Q2n+1' are 
obtained (Step SI 4). 

Then, a curve Ce passing the modified shape-defining points QO 1 , 
Ql 1 , ... , Q2n+l f is obtained based on the position data of the modified 
shape-defining points QO', Ql 1 , ... , Q2n+1' and the first-order differentiated 
values Q0 (1 \ Ql (1> , ... , Q2n+l (1) ' (Step S15). 

For instance, since position data of the two shape-defining points Qj', 
Qj+l f and the first-order differentiated values Qj (1)f , Qj+l (1)t at the respective 
shape-defining points Qj', Qj+1' are given, a third-order curve such as a spline 
curve or a NURBS curve passing these points can be defined. More 
specifically, in the case of defining a third-order curve by a spline curve 
connecting the two points Qj' and Qj+1, coefficients A, B, C and D in the 
following equation representing the spline curve are determined based on the 
position data of the two points Qj' and Qj+1' and the first-order differentiated 
values Qj (1)l and Qj+l (1 \ and thus the three-order curve Ce is created to 
connect the two points Qj ? and Qj+l f (Step S15). 

f^A^ + B^ + Ct + D 

where A, B, C and D are coefficients and V is a curve parameter 
ranging from 0 to 1 .0 

f(t), A, B, C and D are vectors having values with respect to each axis 
of x, y and z. 



t 



r 7 • 

» 

In the foregoing embodiment, an interpolation ratio, i.e. the 
interpolation ratio of the interpolation points (Q2i+1, Q2i+2) between the line 
segment (Pi, Pi+1) is the distance between the line segment (Pi, Pi+1): the line 
segment (Q2i+1, Q2i+2): the line segment ((Pi, P2i+1) = 1: 0.7 : 0.15. The 
interpolation ratio may be other value. 

According to the present invention, a smooth curve approximating an 
original target curve within a tolerance set to the original target curve is 
obtained based on a sequence of command points and interpolation for 
machining is performed on the obtained curve. 



